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In this paper we consider the dynamic Tsallis entropy and employ it for four model systems: 
(i) the motion of Brownian oscillator, (ii) the motion of Brownian oscillator with noise, (iii) the 
. fluctuation of particle density in hydrodynamics limit as well as in (iv) ideal gas. We show that 

the small value of parameter nonextensivity < q < 1 works as non-linear magnifier for small 
values of the entropy. The frequency spectra become more sharp and it is possible to extract useful 
£Nj . information in the case of noise. We show that the ideal gas remains non-Markovian for arbitrary 

^ 1 values of q. 
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I. INTRODUCTION 

There is a great interest towards to studying complex systems of physical, chemical, biological, physiological 
and financial origin by using different approaches and concepts 0,0]. Significant role in these investigations 
plays concept of information en trop y such as Shannon and Renyi entropy, Kolmogorov- Sinai en trop y rate 0, 
S ES 0|- In the Refs. [23, l2lll22i] the notion of dynamical information Shannon entropy |'26| has been 
defined for complex systems. It was suggested to generalize the Shannon entropy by considering square of time 
correlation function (TCF) as a probability of state. Then the entropy becomes a function of time. The dynamic 
Shannon entropy has successfully used to obtain new information in dynamics of i?_R-intervals from human ECG, 
£H , physiological activity of the individuals and short-time human memory |2Ct l2ll |22( • The investigations reveal 
Q 1 also the great role of frequency spectrum of dynamic Shannon ent rop y. 

A new approach for entropy was suggested by Tsallis in Ref. (the earlier works on this subject see in 
Ref. ;2j). This entropy is characterized by one parameter q which was called as nonextensive parameter. The 
I ' natural boundary of this parameter is unit. It means that for q — > 1 the Tsallis entropy transforms to Shannon 
^ . entropy. Nevertheless, it was often considered the cases with q > and even the limit q — ► oo. There is/was 
^\ ' huge activity for application of this new notion to various complex systems (see, for examples, Ref. ^7) and 
, reference therein). It was shown that this entropy gives a new kind of distribution. Many complex systems in 
■ nature are described by this non-Gaussian distribution with different values of parameter q. 

The goal of this paper is generalization of dynamic Shannon entropy in the manner as Tsallis entropy gener- 
• alizes the Shannon entropy. In other words we consider the Tsallis entropy with square of TCF as probability of 
| state. We refer this entropy as dynamic Tsallis entropy (DTE). The real systems in nature have very complicate 
. behavior - the useful signal is lost in "see" of noises and random unexpected influences. For this reason we 
' would like to consider the dynamic Tsallis entropy for model systems, which nevertheless have deep physical 
£h , contents. Application of the present theory to the alive systems was considered in Ref. j2^|. We employ the 
information Tsallis entropy and its frequency spectrum with different values of parameter q for model systems. 
^ \ We exploit four well-known models. First one is the TCF of position of oscillator which performs the Brownian 
motion. In order to consider more real situation we use the motion of Brownian oscillator with noise in second 
model. The noise is modelled by generator of random numbers. Third model is the TCF of relaxation of particle 
density fluctuations in hydrodynamics limit (the Landau-Placzck formula). We consider a specific medium - 
Helium. The fourth model describes the relaxation of particle density fluctuations in ideal gas. In the later 
case we additionally calculate the spectrum of non-Markovity parameter for different q. The non-Markovity 
5^ ' parameter and its spectrum firstly was suggested in Refs. |14l Il5l| and it was calculated for ideal gas in Ref. 
[l9|. This parameter characterizes the statistical memory effects in systems. In present paper we use slightly 
different approach and define the non-Markovity parameter in terms of the entropy. 

In all considered models the behavior of the dynamic Tsallis entropies and their frequency spectrum crucially 
depend on the parameter q. We observe that decreasing of this parameter q — * allows us to enlarge the fine 
structure of entropies and to make their frequency spectrum more sharp. 

The organization of this paper is as follows. In section [H] we discuss the well known hierarchy of Zwanzig- 
Mori's kinetic equations and general properties of the DTE and define different kind of relaxation times. The 
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TCF of Brownian motion of oscillator is considered in section II I II and with noise in section IIVI The formula 
of Landau and Placzek for relaxation of density fluctuations are exploited in section [V] The TCF of density 
fluctuations in ideal gas is used in section IVTl We finish our paper by concluding remarks in last section IVTT1 



II. THE ZWANZIG-MORI HIERARCHY AND DYNAMIC TSALLIS ENTROPIES 

At the beginning we shortly discuss the well known hierarchy of Zwanzig-Mori equations 0, |24| . Let us 
consider the dynamic variable 6A(t). It may be, for example, the Fourier component of density fluctuation. 
This variable obeys to the Liouville equation: 

= iL5A(t). 

Applying n times the Liouville operator L to the initial variable SA(0) we obtain the infinite set of variables 
B n (p) 

B n (0) = L n 6A(0), 

by using which and the Liouville equation we may obtain the initial dynamical variable in arbitrary moment of 
time 



m 

n—Q 

Applying the Gram-Schmidt orthogonalization procedure to this set of functions we obtain the complete 
set of dynamic variables W n which are orthogonal at the initial time 

(w^Wt ) = (\w n \ 2 )s 7hh 

where (. . .) denotes the statistical average over Gibbs ensemble. If the dynamic variable is evaluated by the 
Liouville's operator, then the orthogonality is preserved at any moment of time due to the sclf-adjointness of 
the Liouville's operator. 

The time correlation function Mq of variable Wq = 8 A is defined as follows: 



M (t) 



(W (0)*Wo(t)) (Wo* exp (iLt)W ) 



(\Wo(0W) (|Wo| 2 ) 
It is well known that this TCF obeys to the infinite hierarchy of the Zwanzig-Mori's kinetic equations 

= i4 n) M n (t) - n 2 n+1 j* dt'M n+1 (t')M n (t t% 

where 

i(n) _ {W*LW n ) Q2 = (\W n \*) 



(1) 



ui, 



° (l^| 2 ) ' n (|W„_n 
and 



_ (W*exp(dj n 2 k)W n ) 



The operator L 1 ^ is defined in following manner: 

L$ = Pn-lPn-2 • ■ • PoLPq . . . P n - 2 P n -l 

in terms of projectors P n = 1 — II n , where the II n is projector for state W n : 

_ w n )(w: 



3 



The functions M n for n > 1 are not, in fact, usual TCFs because the operator exp (iL^t) is not operator of 
evolution. Nevertheless, we will refer for these functions as TCF's for the next dynamic variables. 

The functions M n are considered as functions characterizing the statistical memory of the system. In order 
to describe quantitatively the non-Markovity of hierarchy, the parameter of non-Markovity and its spectrum 
were introduced in Refs. [bH Il5| . The spectrum of this parameter is defined as ratio of two neighboring 
relaxation times 

e n = —■ (3) 



The relaxation time r„ of function M n was defined as real part of the Laplace image (see Eq.(0}) of this function 
at zero point: 



3?M„(0) = 5? / M„(t)dt 



o 



where 3? means a real part. Because of the function M n+1 is an integral core of integro-differential equation 
JPl for M n , then this parameter © compares the integral core with function itself. More precisely we compare 
squares under TCFs. If parameter e„ is around unit, this level (level means n) is non-Markovian: we can not 
transform integro-differential equation for M n to differential one, and vice versa, if this parameter tends to 
infinity the core has sharp peak and we may transform the integro-differential equation for M n to differential 
one. In this case there is no memory (integral) in this level. In this paper we use slightly different definition for 
relaxation time (see Sec lVIf) . 

Applying the Laplace transformation 



M n (s)= / dte- st M n (t) (4) 
Jo 

to equation Q we transform this hierarchy to the infinite system of algebraic equations 

M n (s) = [s- *w< n > + n n+1 M n+1 {s)}-\ (5) 

by using which we may express M n in term of zero TCF Mq. 

In statistical physics of non-equilibrium systems the time correlation function (TCF) acts as the function of 
distribution and pair correlation and can be used to calculate different thermodynamical parameters and the 
spacial structure of the system . For many physical discrete systems it is impossible to find a distribution 
function. For this reason it seems optimal to obtain the TCF for investigations of complex systems with the 
help of integro-differential equations which are based on small increments of time and independent variables. 

The set of the measured parameters of the complex system may be represented as a set of fluctuations [2^ 

Z = {C(T), C(T + r), C(T + 2r), ■ • • , £(T + fcr), • • • , C(T + tN - r)}, (6) 

where £ represents the fluctuation Sx of some quantity x. 

Let us define the sampling by length k, which starts at the moment T + mr by the relation 

Of* = K(T + mr), ■ ■ ■ , C(T + mr + (k - I)r)}. 

The operator which projects the sampling by length k on the sampling at the initial moment of time has the 
form 



n = 



where angle brackets mean the scalar product (time average). The utilization of this projection operator allows 
us to represent the sampling as a sum of two independent parts 



where 



C+k = n C+fc = C° k M (t), (7a) 
C+fc" = (i - H)C+* = C+fe - CkMo(t), (7b) 



((Cm+fe) ) — ((Cm+fe') ) + ((Cm+fc") )) 



and TCF M Q (t) is defined by ©. 
It is easy to see that 



which is the consequence of orthogonality of 17a|) and (|7b|l . 
Direct calculations yield 

((c +fc o 2 > = m) 2 )M (tf. 

For stationarity processes, when dispersion does not depend on time, we obtain 

((C +fe ") 2 > = <(C +fc ) 2 > - ((a°) 2 )Mo(i) 2 - ((C +fe ) 2 )(l - Afo(t) 2 ). 
Therefore the mean-square value of the fluctuations is presented as a sum of two parts 

<(C +fc ) 2 > = <(4°) 2 )AfoW 2 + <(C£) 2 )(1 - AfoW 2 ), (8) 

or in a generalized form 

<(C +fc ) 2 ) = <(4°) 2 )Mo« 2 + <(C +fc ) 2 >(l - Mo(tf)- (9) 
The above expression @ has a standard form of a mean value 

<(C +fc ) 2 > = ((4°)>(*) + <(C +fc ) 2 )QW, (io) 

where P(i) is the probability of the state and Q(t) = 1 — -P(t). By analogy with this formula we define 

P n (t) = \M n (t)\ 2 (Q n (t) = 1 - |Af„(t)| 2 ), n > (11) 

as the probability of the creation (annihilation) of correlation of fluctuations (or memory) for the nth level of 
relaxation (see, for details Refs. uM) . 

In the recent work of authors |22j the dynamical informational Shannon entropy for the study of complex 
systems was suggested: 

S n (t) = - PWln J P,(t) = -|Af„(t)| 2 ln|M„(^)| 2 -{l-|M n (t)| 2 }ln{l-|Af„(0| 2 }, (12a) 

i=c,a 

where 

S c n (t) = -\M n (t)\ 2 \n\M n (t)\ 2 , (12b) 

S a n {t) = -{l-\M n {t)\ 2 }\n{l-\M n {t)\ 2 }. (12c) 

Here the S^(t) is the entropy for the stochastic channels of memory creation, and S%(t) is entropy for the 
stochastic channels of memory annihilation. 

The infinite set of TCFs M n (t) produces the infinite set of entropies which are defined by relations 

S n [t] = S[M n (t)], S nc [t] = S c [M n (t)}, S na [t] = S a [M n (t)}- 

Because of the TCFs are always smaller then unit then all of these entropies are positive. 

We also generalize our definitions of entropy due to Tsallis by introducing parameter of nonextensivity q: 

S q (t) = J^\^-l + \Mo(t)\ 2 ^ (13a) 



sm = JW)\*-\W)\\ 

sa{t) = Jl-\Al(m«-l + \M (t)^ 
q-l 



The infinite set entropies are defined by relations 

S?[t] = S q [M n (t)}, S?[t] = S c q [M n (t)}, S™[t] = S a q [M n {t)}. 
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We also calculate the frequency spectra S n [v] , S nc [v] , S na [u] , S q [u] , S q c [v] , S q a [u] of these functions and define 
the spectra of the relaxation times by relations 

r n = S n [i>] v =0) T nc — <5" lc H^=o, T na — S na [v] v=0 , 

Tqn = Sg[y] u —0, T q71 c = Sq C [v] v= Q, T qna = S q a [v] v= Q. 

We define the frequency spectrum as Fourier transformation of these quantities by relation 

r+oo 



M n {v) = / M n {r)e-^ WT dr. 

J — oo 

To calculate the Fourier transformation of all quantities we use package for numerical Fourier transformation 
from package " Mathematica" . We tabulate functions M n (r) with step A = 1/40 in range r = (0,60). The 
amount of points N = 2401. Then the Fourier transform 

/+oo r + co 

M„(T)e^™ T dT = 2 / M„(t) cos(27WT)dT 
-oo 

N-1 







2A5R J2 M n {kA)e~ 2 * lkv/N - M„(0) 



k=0 

is represented as a function of discrete frequencies v — j^g (k = 0, . . . , N — 1) 

N-1 

M n {v) = 2A3?^ M n (kA)e- 27rUk / N ~ M n (0). (14) 
;=o 

Let us summarize here some analytical results. In the limit q — > 1 the formulas Ijl3(l transform to the equations 
l)12[l. For this reason we may refer for entropies 1)1 2[l as entropies at the point q = 1. The entropy S q Ijl3(l 
amounts to maximum value, (2 1 ~ g — 1)/(1 — q), at the point |M n (t)| 2 = |. For small values of |M„(i)| 2 the 
dependence of entropies S q , S q ° changes drastically on value of g. For g <C 1 we have 

S q « |Afo(0| 29 , 5 9 c « |M (i)| 29 , ^ « |MoW| 2 . (15) 

These expression are valid for all cases: |Mo(i)| 2 <C q,\Mo(t)\ 2 ^> q and |Mo(t)| 2 ~ q. For comparison we 
have for small values of \M {t)\ 2 : S w -|M (i)| 2 ln|M (t)| 2 . Therefore < S°,S q ,S and we observe that 
decreasing of <7 does as magnifier for small values of \Al (t)\ 2 . 

For q » 1 and still |M (i)| 2 « 1 we have to consider three cases q\M (t)\ 2 <C 1, q\M (t)\ 2 ~ 1 and 
q\M (t)\ 2 ^S> 1. We have correspondingly 

S g « |M (t)| 2 , S g c «i|M (t)| 2 , S««|M (i)| 2 , [g|M (i)| 2 « l] 

5 g « _±((i_|M (i)| a )«-l), S g c « Vo(t)| 2 , S««-i((l-|M (f)| 2 r-l), [<?IW)| 2 ~ l] 
q q q 

S q « -, S| « i|M (i)| 2 , [ 5 |M (t)| 2 » 1] . 

q H q H q 

We observe that for sufficiently large g and small |M (i)| 2 Clwe have: S q , S q , S q <C 5*. Therefore we observe 
that increasing of q do as demagnification lens for small values of |Mo(t)|. We note that the entropies equal to 
zero for arbitrary q for zero value of M (t) as well as for M (t) = 1. It is easy to see, if the TCF M (t) possesses 
an extremum at some time to the entropies have extremum, too. The structure of extremums of TCF generates 
the same structure of extremums of entropies. After this general frameworks let us consider specific models. 



III. THE MOTION OF BROWNIAN OSCILLATOR 



Let us consider the Brownian oscillator which models the particle with internal oscillatory degree of freedom 
with frequency loq. The friction coefficient 8 describes the movement of the particle in a medium. We may use 
the general theory of Zwanzug and Mori [Tol |24| for this particle because its position and velocity are random 
quantities. 
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The position and velocity of the Brownian oscillator is subject for equations 

I = «' ^ 
dv a , F(t) + f(t) 

- = -2cv-u> x + , (17) 

where c = (3 /2m, -F(i) - external force, /(t) - the Langevin random force. The random quantities x,v,f are 
described by the following average values: 

(x) = 0, (v) = 0, (/) = 0, 

(.t 2 ) = -^2, (« 2 ) = -, (xv) = 0, 
mujQ m 

(xf) = 0, (vf) = 0, (f(t)f(t'))=2Tl36(t-t'), 

where T - temperature in the units of the energy, m - the mass of Brownian particle. 
We use the fluctuation of particle's position x 

Wo = x(0) 

as the initial dynamic variable. The next orthogonal dynamic variables are calculated with help of the following 
recurrent relations: 



Wi = UL-u, 
W„. = 



{iL-4 0) }w Q , 

[iL-w^-^Wn-x+^l^Wn^ n>l, (18) 



where the frequencies Wq™' ) and f2„ are defined by equations 

r„) _ (W*iLW n ) o2 _ (\W n ^ 



{\W n \>) ' " {\W n - 
The straightforward calculations give 

(0) _ (x(0)v(0)) _ 

" (W0)l 2 > _u ' 

Wi = (iL-lu ( 0) )Wo =v(0), 



W ~ /IW |2\ ' n n~—r ia\- ( 19 ) 



w 2 



(H0)\ 2 ) 

(\V(0)\ 2 ) 2 

(|x(0)P) °' 
/(0) 



-2c, 



m 

Taking into account the equations of motion {T7J|, one obtains the time correlation function (TCF) 

(afoMO)) r-e^ - r+e^ 

Mo{t) ~ TRoFT " — — ' (20) 



where r± = — 



To find the TCF for next levels we exploit the Zwanzig-Mori infinite chain of equations given by Eq. 
The Laplace-image of initial TCF has the following form: 

Mo(s)~ -< r - +r +> 



(s - r_)(s - r+) 

With n = we obtain from chain |£| the Laplace-image of first memory function 

Mi{s) = 1 ■ 

s — (r_ + r+j 
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Taking the inverse Laplace transform we obtain the first memory function 

Mi(t) = e ( r - +r +)* = e~ 2ct . 

The calculation in closed form of the memory functions on the next relaxation levels requires the following 
detailed elaboration of the structure of Langevin force /. For example, in order to calculate the memory 
function M2 (t) we need for frequency f^2 which may be found by setting the mean value of square of Langevin 
force. For this reason we restrict our consideration by first two levels of relaxation. 

We consider the TCF of Brownian oscillator in the case of small damping p — c/ujq « 1. In this case the 
TCF M (t) has the following form: 

M «< ! '=smT =cos(wt)e " 01 " <2i) 

and describes the motion of Brownian oscillator with frequency u> = 2ttv' ^> c = (3 /2m, where m is the mass of 
oscillator, and (3 is the friction coefficient of Brownian particle [l2|. 

There are two parameters v 1 and c which characterize the frequency oscillation and relaxation damping of 
TCF, respectively. For simplicity we consider dimensionless time t, frequency of oscillation v 1 and damping 
parameter c. The Fourier transform of square of this TCF has the following form 

— — — n r> r> 

f2 



M 2 {v) 



2(c 2 +tt 2 v 2 ) 4(c 2 + ir 2 {v -2v'Y) A{c 2 + tt 2 {v + 2v') 2 ) 



We consider frequency spectrum of square of TCF because the entropies are expressed in terms of square of 
TCF. Therefore there are three maximums at v — 0, v = ±2z/. In the limit of zero damping c — > we obtain 



by using well-known formula: 



M » = ~8{v) + \8{v 2v') + ls(u + 2v') 



tt5(x) = lim 



ar 



In this case the spectrum consists of three lines at v = 0, v = ±2v' . 

In Fig we reproduce the entropies l|12l) and in Fig|2] we show the spectra of them for v' = 0.1 and for 
c = 0, 0.01, 0.1, 1. There is a peak in spectra at double frequency v = 2v' = 0.2 for arbitrary small but non-zero 
damping c. For zero damping c = this peak disappears in total entropy. There are peaks in S c and S a at 
v = 0.2 but with opposite sign. We note that the increasing of the damping leads to smearing the fine structure 
of entropies. 

Calculation the entropies l|13fl shows that the small values of q "works" as non- linear magnifier for small value 
of Mq . To illustrate this fact we reproduce in Fig. [3]the entropies (|13fl and (|12fl for 5=1, 0.1, v' — 0.1, c = 0.1. 
First of all, as expected from Eq. I)15|l the smaller entropy the greater magnification. Second, the entropy 
does not change sufficiently. The great variation prove S q and S^. The decreasing of q leads to increasing 
the value of peaks for great frequencies and makes peaks more sharp. The small q makes better sharpness of 
frequency spectra. This is not the case for 5°. It does not change great. 



IV. THE MOTION OF BROWNIAN OSCILLATOR WITH NOISE 

The real signal from alive systems often contains noise (see Ref. |2^|). For this reason we suggest the model 
of Brownian oscillator with noise. We consider the following model of TCF: 

M (t) = R(t) cos(27iVi)e~ c|t| , (22) 

where R(t) denotes the random numbers in interval (—1, 1) and R(0) — 1. Therefore the function R(t) makes 
random the amplitude of oscillation but this is not the case for frequency and damping parameters. At the 
beginning we know the frequency of oscillations. The TCF has more complicate form which is more close to real 
dates. The random numbers describe a noise which usually appears in an experiment. The time dependence of 
entropy is much more complicate but nevertheless in this case the DTE works better. We observe the appearing 
peaks in places, which we know the peaks must be, but they disappeared in noise. There is another observation: 
the noise is better for frequency spectrum. The peaks for small frequencies looks better (see Fig. 
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FIG. 1: The plots of entropies (II 12 P for frequency v' = 0.1 and c = 0(a), 0.01(6), 0.1(c), 1(d) [y' and c are oscillation and 
damping parameters of TCF 1211 1. It is seen that the time behavior of total S and single channels S c and S a of DTEs 
for various damping regimes reveal a stochastic ordering of time correlation. 



DENSITY FLUCTUATIONS IN HYDRO DYNAMICAL LIMIT 



The TCF of density fluctuations in hydrodynamical limit was calculated by Landau and Placzek [l^. It 
describes scattering of light in liquid in hydrodynamics limit when k — > 0. The TCF has the following form 



M (t) 



cosd s kt + (1 — a)e 



-a'k 2 t 



where 7' = ^ ^77 + ( + re 



a = c v /c p . Here the c v ,c p , k,7],£, p,"& a are specific heat 
capacities in units of mass at constant volume and constant pressure, the coefficient of thermal conductivity, 
the coefficient of shear viscosity and volume viscosity, the mass density and sound velocity, correspondingly. 

The spectrum of this TCF contains three peaks. The central Rayleigh peak at zero frequency describes 
isothermal propagation of sound. Two symmetric peaks at frequencies to = ±z9 s fc describe adiabatic propagation 
of sound with damping (Brillouin doublet). 

It is more suitable to define dimensionless time r by relation r = , d s kt/2'K. Then the position of Brillouin 
doublet will be at the dimensionless frequency v — uj/2tt = 1, and TCF will take the following form 



M (r) 



cos 2-7TT + (1 — ct)e 



— ok2-KT 



where 7 = 2^?: (l 7 ? + C 



a 



c v /cp. We consider a specific medium - Helium 



at temperature T = 20°C and pressure p 



6 ■ 10 9 m 1 , (7 : 



pc p S s ' 

16. In this case we have ^8|: a ~ 0.56, i9 s ~ 272 m/c, 7 



7 • 10 m . We make calculations for k = 2 • 10 m 




-0.3 L -0.3 



FIG. 2: The plots of spectra of entropies 11211 for v' = 0.1 and c = 0(a), 0.01(&), 0.1(c), 1(d). One can see that a stochastic 
ordering of time correlation in previous FigQ reduces to an appearance of specific peculiarities in low frequency region. 
Weak damping (c = 0.1) results in amplification of characteristic frequency peaks at v — 0.2 and v — 0.4 whereas zero 
damping (c = 0) leads to a dissapearance of specific frequency peaks. 

In the Fig[S]we reproduce time and frequency dependencies entropies defined before for two values of parameter 
q = 1,0.1. There is periodicity over r with unit period which gives the appearance specific peaks in the frequency 
spectrum (Brillouin doublet). We observe the same picture as in previous section. Decreasing the parameter of 
nonextensivity q leads to increasing small peaks in entropies S qn and S qcrrin , whereas quantity S qamn is changed 
insufficiently. The frequency spectrum entropies S qn and S qcrrin becomes more sharp. The shape of frequency 
spectrum S qa m n , in fact, does not change. 

The non-Markovity parameter of this system has been calculated earlier in Refs. It was shown that 

in hydrodynamical limit k — > the spectrum of non-Markovity parameter has a form of alternating Markovian 
and non-Markovian levels. 



VI. IDEAL GAS 

Let us consider the Fourier transformation of the fluctuation of the particle number density of a system 

1 N N 
$Pk(t) = — exp(ikri) - tf4,o- 
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FIG. 3: The plots of entropies (1131 for q = 0.1 and 11211 (q — 1) and their spectra for v' — 0.1, c = 0.1. Comparison 
of dynamic Shannon (q — 1) and Tsallis (q — 0.1) entropies shows the amplification of DTE for small values of TCF. 
Therefore one can conclude that DTE acts as a magnifier for small values of TCF. 
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FIG. 4: The plots of entropies il'-ii for q = 0.1 and il l 21 (q = 1) and their spectra for v' = 0.1, c = 0.1 and for Brownian 
oscillator with noise. Comparison of dynamic Shannon (q = 1) and Tsallis (q = 0.1) entropies shows the amplification 
of DTE for small values of TCF. One can conclude that DTE acts as a magnifier for small values of TCF. 



For this case the initial TCF is calculated exactly and it has the following form 12]: 

(Spk(t)6pk(Q)*) =e -t 2 /t 2 r f 2 = n ^/,„2r 



M (t) = 



(\S Pk (0)\ 2 ) 



, tf. = 2m/k z T. 



For this variable the all frequencies Uq are equal to zero. The main relaxation frequencies in Eq. Q have 
the simple form 



Q„ = ntt{, tt{ = k 2 T/m. 



11 




Let us rescale time t — > r = tS7i , and the Laplace transformation parameter s — > c = s/Qi, and the Laplace 
images M n (s) — > fh n (s) — QiM n (s). In this case the hierarchy (|TJ has the following form 



»« n+ i(c) = 



1 



1 



n+1 [ m n {c) 
The Laplace image mo(c) may be found in close form: 



mo(c) = e 2 



(23) 



where Erfc(x) = 1 — Erf(a:) - additional probability integral. 

The inverse transformation may be represented in the following form: 



m «(t) = — / e CT m n (c)dc, 

■^ 7r * J a—ioo 



where c is greater then real part of zeros of m n (c). By using the expression (|23|l we may set cr = and by 
changing c — * ix we obtain 



1 /"+ 00 

M n( T ) = 7^ / e icx m n {ix)dx, 



or in manifest form 



Mi(r) = 



2tt 



1 



1 / >+0 ° 



mo(ix) 
1 



dx, 



dx. 



These formulas we will calculate numerically. The first three memory functions and their spectrum are plotted 
in Fig El In Fig. 0we reproduce the plot of the entropy S n and their frequency spectrum for n = 0, 1, 2. 

It is not difficult to show that the entropy S n (|12|l amounts to maximum value, In 2, at the point |M„| 2 = I. 
For n = the position of maximum is at the point r = Vln 2 « 0.832. The greater n, the smaller time of 
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FIG. 6: The first normalized memory functions M„(t) and their frequency spectrum M n (y) for ideal gas. 





FIG. 7: The plot the time dependence of entropy and its frequency spectrum for different relaxation levels n — 0, 1, 2 for 
ideal gas. 

maximum. The relaxation times To ~ 1.6453, t\ — 1.0688, T2 = 0.8630. The parameter of non-Markovity, 

0E3 

has the following values eo ~ 1.54, ei ~ 1.24. Note that these values very close to that calculated in paper [l^ | 
directly for TCFs, e ~ 1.57, e x ps 1.27. 

In the Fig. JSJ we give the plots of entropy S™ and their spectrums for q = 1/2 and q = 3/2. In order to show 
what happens if we will vary the value of q we reproduce in Fig. the entropy and its spectrum for n = 2 
and for q = 0.5, 1, 1.5. To show the dependence for more wide range of q in Fig. JTJJJ we give entropies for 
more wide range of q. 

By using these three quantities we define seven different spectra of non-Markovity parameter: 



' qan 



Tqcn 
7~qan-\-l 



Tqcn-\- 1 
T, 



1 C <?O.Tl 



") €qac n 



Tqan-\-l 

qan T~a 
) €qcn 



1 qcn 



' qcn+1 



Tqn+l 



i tqan 



1 qan 
Tqn+1 



(24a) 
(24b) 



The plots of all quantities are shown in Fig. Qllfl. As expected at the beginning more interesting situations 
are possible for small values of q. For great value of q all lines tend to constants. 



VII. CONCLUSION 



Let us summarize our results. We considered the dynamic Tsallis entropy and applied it for model systems. In 
accordance of Refs. [13, mm the square of TCF is regarded as probability of dynamic state. We used the time 
dependent entropies as well as their frequency spectra. We considered four model TCFs which, nevertheless, 
have physical sense. The first and second models describe the motion of Brownian oscillator without and with 
noise. The noise is modelled by generator of random numbers. Third model is the Landau-Placzek model 
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FIG. 8: The plots of time dependent entropy Sq and their frequency spectrums for q = 1/2 and q — 3/2. 
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FIG. 9: The time dependent entropy Sq and their frequency spectrum for n = 2 and for q — 0.5, 1, 1.5. 

for TCF of particle density fluctuations in hydro dynamic limit. The last model is ideal gas, the relaxation of 
particle density fluctuations. 

In all model we considered different values of parameter of nonextensivity q. All entropies have the same 
structure zeros and extremums as square of TCF. The magnitude of extremums sufficiently depends on value of 
q. Small values of q work as non-linear magnifier: the smaller magnitude of entropy the greater magnification. 
Great values of q do in opposite way as demagnified lens. Concerning the frequency spectra of entropies we 
observe that for small values of q the peaks become more sharp and more large. For systems with noise this 
property works better. It is possible to reveal peaks even if they are lost in noise. 

For ideal gas we additionally calculated the spectrum of parameter non-Markovity by using different definition 
of relaxation time. By using three kind of information Tsallis entropy we defined seven kind of spectrum of 
parameter non-Markovity. We observe that all of these parameters (except e gca „) are close to unit. It means 
that the ideal gas remains non-Markovian system for arbitrary value of q. It is in qualitative agreement with 
Ref. 0. The variation q from unit does not mean the appearance of new interactions. The system is became 
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FIG. 10: The time dependence of the entropy Sq and its frequency spectrum for six values q — 0.001, 0.1, 0.5, 1, 1.5, 10 
and for n = 0. One observe that the increasing Tsallis' parameter q leads to trample a quantity of S q in domain of short 
time. Due to this fact one can make sufficient amplification in domain of short time owing to variation parameter q. 

ideal and non-Markovian from this point of view. 

Our analysis allows us to conclude, that use of dynamic Tsallis's entropy extends essentially possibilities of 
the stochastic description of model physical systems. Advantage of application of DTE is that it allows to 
strengthen or suppress fluctuations either in low-frequency, or in high-frequency areas of a spectrum. Similar 
supervision opens appreciable prospects in the field of the study of real complex systems of wildlife where 
dynamic states of physiological and pathological systems are very important. 
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FIG. 11: The spectrum of all parameters of non-Markovity as function of q. The influence of parameter q is most 
effectively for small values of q <C 1. In all cases we observe strong non-Markovity and clearly marked effect of statistical 
memory. The only case where we may sufficiently increase the non-Markovity properties by decreasing parameter q — > 
is cross term ca (see Eq. <!24bll ). For small enough q <C 1 the non-Markovity parameter e qcarl may reach great value 
e gc a„ > 1. 
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